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ABSTRACT 

It is demonstrated here that if the prompt GRB emission is produced by the 
simplest version of the external shock model, a specific relation should prevail 
between the observed duration, isotropic equivalent energy, and photon peak en- 
ergy. In essence, this relation arises because both the burst duration and the 
typical energy of the emitted synchrotron photons depend on the same combi- 
nation of the, usually poorly constrained, external density at the deceleration 
radius, n^ec, and initial bulk Lorentz factor, r . This has the fortunate con- 
sequence of making the relation independent of both T and n^ec- Unless the 
efficiency of electron acceleration is very low, synchrotron gamma-rays from the 
external shock would fail to meet the current observational constraints for the 
vast majority of GRBs, including those with a smooth, single peak temporal pro- 
file. This argues either against an external shock origin for the prompt emission 
in GRBs or for changes in our understanding of the microphysical and radiation 
processes occurring within the shocked region. 

Subject headings: gamma-rays: bursts - hydrodynamics - ISM: jets and outflows 



1. Introduction 

The simplest version of the standard fireball model for gamma-ray bursts (GRBs) in- 
volves a spherical explosion taking place in a uniform or a stratified surrounding medium. 
When an explosion deposits a large amount of energy into material with a much smaller 
amount of rest energy within a compact volume, an ultra-relativistic pair fireball is formed 
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(Cavallo & Rees 1978; Paczynski 1986; Goodman 1986). The large pressure of the explosion 
causes the fireball to expand, and the thermal energy of the explosion is transformed into 
bulk kinetic energy due to strong adiabatic cooling of the particles in the comoving frame. 
Because of the Thomson coupling between the particles and photons, most of the original 
explosion energy is eventually carried by the baryons that were originally mixed into the 
explosion (Shemi & Piran 1990). This bulk kinetic energy cannot be efficiently radiated as 
gamma rays unless it is converted back to internal energy (i.e. the velocities of the protons 
must be re-randomized). This requires shocks, and in order to tap a reasonable fraction of 
the total kinetic energy, the shocks must be (at least mildly) relativistic. 

Impact on an external medium would randomize about half of the initial energy merely 
by reducing the expansion Lorentz factor by a factor of ~ 2 (Rees & Meszaros 1992). Al- 
ternatively, internal shocks may form within the outflow: for instance, if the Lorentz factor 
of the outflow varied by a factor > 2, then the shocks that developed when fast material 
overtakes slower material would be internally (at least mildly) relativistic. There is a general 
consensus that the longer complex bursts must involve internal shocks, while simple smooth 
profiles could arise from an external shock interaction (Sari & Piran 1997a,b; Ramirez-Ruiz 
& Fenimore 2000; Nakar & Piran 2002; McMahon et al. 2004; Ramirez-Ruiz & Merloni 2001). 
The latter would in effect be the beginning of the afterglow. 

An external shock moving into a medium with a smooth density profile would naturally 
result in a burst with a simple time-profile. Angular variations within the outflow might still 
cause variability in the light curve, but variations on very small angular scales (9 < T^ 1 , 
where Tq in the initial Lorentz factor) are required in order to produce the large variability 
of the prompt GRB emission (Fenimore et al. 1999; Dermer & Mitman 1999). A blobby 
external medium could produce significant variability only if the covering factor of blobs 
is low, implying modest efficiency Furthermore, the resulting variability in the light curve 
would be small if produced close to or after the deceleration radius, or if the portion of the 
ejecta that collides with a blob is decelerated significantly (Nakar & Granot 2006). 

The purpose of this Letter is to demonstrate that if the prompt emission is produced 
by the simplest version of the external shock model, this implies a specific relation between 
the observed duration, isotropic equivalent energy (or luminosity), and photon peak energy, 
which is apparently incompatible with observations. This relation is derived in § 2 and 
compared to observations in § 3. The implications are discussed in § 4 along with possible 
caveats. 
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2. External Shock Model 

In the simplest version of the external shock model, the outflow is approximated by a 
uniform thin shell. A forward shock is driven into the external medium by the outflowing 
ejecta, while the latter is decelerated by a reverse shock (and/or by pdV work across the 
contact discontinuity that separates it from the shocked external medium). The dynamics of 
a spherical shock wave eventually approaches a self-similar evolution (Blandford & McKee 
1976) which depends only on the explosion energy E and on the external mass density 
Pcxt = n cx tm p (the Lorentz factor depends only on their ratio, E/p ext ). If the initial GRB 
outflow is collimated, an additional parameter - the jet initial half-opening angle, 9 , is 
required in order to specify the flow. However, for r #o 3> 1 (as appears to be the case from 
afterglow modeling; Panaitescu & Kumar 2002) the dynamics at early times - before the jet 
break time (as long as T > Oq 1 ) do not significantly deviate from the spherical case, where 
the true kinetic energy E is replaced by its isotropic equivalent value E iso . Therefore, it is 
still valid to adopt the spherical dynamics for the prompt emission from the external shock 
in this case as well. 

Most of the energy is transfered to the shocked external medium at the deceleration ra- 
dius, -Rdcc, where the inertia of the swept-up external matter starts to produce an appreciable 
slowing down of the ejecta. For a given shock dynamics, the luminosity and spectrum of the 
emitted radiation are determined by the fractions e B and e e of internal energy in the shocked 
fluid that are carried, respectively, by the magnetic field and by relativistic electrons, as well 
as by the shape of the electron distribution function. 

As seen in the rest frame of the downstream fluid, most of the mass and of the kinetic 
energy of the incoming upstream fluid is in protons (or other ions), unless the external 
medium is highly enriched in pairs. Therefore, a simple isotropization of the velocities 
of the upstream particles at the shock transition would give the electrons only a very small 
fraction of the total internal energy (~ m e /m p ). This would imply a very small radiative 
efficiency, since the radiation is emitted primarily by electrons. For a radiatively efficient 
system, physical processes must therefore transfer a large fraction of the swept-up energy 
to the electron component. The energy of the particles can be further boosted by diffusive 
shock acceleration as they scatter repeatedly across the shock interface, acquiring a power 
law distribution dN e /d^ e oc % p at 7 e > j m . 

The strength of the magnetic field is another major uncertainty. Most of the required 
magnetic field must typically be generated in-situ, presumably through plasma instabilities or 
turbulent motions, but its strength has yet to be derived from first principles. The standard 
prescription is to assume that the magnetic field energy density U' B = (B') 2 /8ir is a fixed 
fraction e B of the downstream proper internal energy density, B' = (32ireBn ext m p c 2 T 2 ) 1 / 2 , 
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where primed quantities are measured in the comoving frame. 
The typical (minimal) electron energy is given by 

m p I p — 2 \ r, 



7T7 



m e \p - 1 J £ e 



(1) 



where Y is the Lorentz factor of the fluid behind the forward shock, and £ e is the number of 
relativistic electrons (or positrons) per proton, which for a proton-electron plasma is equal 
to the fraction of the electrons that are accelerated to relativistic energies. 1 



The peak synchrotron photon energy is given by 
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where T 2 = r(R^ ec )/100, n is nd ec = n cxt (Rdcc) in units of cm 3 , and g = 6(p — 2)/(p — 1) 
(where g — 1 for p — 2.2). For p cxt = n cxt m p = Ar~ k (with k < 3) we have 
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where 2 r = r(i?d cc ) = Tdec, a = 2a 2 with a<i ~ 1, and 
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In the external shock model the duration of the GRB is T G rb ~ ^dec, and therefore 
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1 It is assumed here that all the relativistic electrons take part in the power law distribution of energies; 
the definitions of e e and £ e would not include possible additional components, such as a thermal component. 

2 Note that the initial Lorentz factor of the outflow can be higher than T (if the reverse shock is rela- 
tivistic). 
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This relatively simple relation between different observable quantities arises since Td ec oc 
R dec /T 2 dcc oc ^(naeerL)- 173 while E p oc TB'-fa oc ef \e e j ' ie) 2 {n Acc Y\ cc ) x l 2 , so that both E p 
and T G rb ~ Tdec depend on nd ec and r dcc only through the combination n dcc rj cc . Therefore, 
the dependence on both nd cc and Tdec (which are hard to determine from observations) can 
be eliminated by taking the combination E P T^ oc e 1 ^ 2 (e e /$, e ) 2 ElJ 2 . 

The strength of Eq. 7 is that it depends mainly on quantities that can either be directly 
measured, like the peak photon energy (E p ) and the duration of the GRB (Tqrb), or that 
can be reasonably constrained by observations. Here -E^o is the isotropic equivalent kinetic 
energy of the outflow, which for a reasonable radiative efficiency, e 7 ~ 0.5, is of the order of 
the isotropic equivalent energy output in gamma-rays, E ltiso , that is measured directly. 

3. Comparison to Observations 

In order to compare the limit imposed by Eq. 7 with observations, we use the following 
observational properties derived by Ghirlanda et al. (2004) and Kaneko et al. (2006): T90, E p 
and .E^iso- In Fig. 1 we show the distribution of ^ as a function of T90, E p and -E-y^so- Filled 
circles are typical long bursts from the sample compiled by Ghirlanda et al. (2004), while the 
empty circles are the four long GRBs found so far to be spectroscopically associated with 
type Ic supernovae (Kaneko et al. 2006). Of the latter, three have a smooth, single peak 
temporal profile (while GRB 030329 has two peaks). Only two bursts have \P < 1 while most 
bursts (and in particular those associated spectroscopically with supernovae) have ^ 3> 1. 
Fig. 2 shows the maximal value of £ e that is allowed according to Eq. 7, £ e ,max = ^~ l l 2 - 

There are some necessary limitations to our approach. The choice ofe e = e£ = l/3 
that has been used in the definition of ^ in Eq. 7 is conservative. More typical values that 
are inferred from afterglow modeling (e e ~ 0.1 and €3 ~ 0.01) would result in the values 
of £ e being smaller by a factor of 8.0(e e /0.1)~ 1 (e^/0.01)~ 1 / 4 when compared to £ eimax - It is 
also important to note that -E 7i i SC , is used as an estimate for the isotropic equivalent kinetic 
energy E iso . This would increase the value of £ eimax by a factor of [(1 — e 7 )/e 7 ] 1//4 , where e 7 
is the 7-ray efficiency: E- lso /E 1 ^ so = (1 — e 7 )/e 7 . However, even for 7-ray efficiencies as low 
as e 7 ~ 10 -2 , £ ei max would only increase by a factor of ~ 3. 

4. Discussion 

It has been shown that the simplest version of the external shock model implies a 
relation between different observed quantities of the GRB (Eq. 7), which can conveniently be 
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expressed in the form £ e < £ eimax = X I / ^ 1//2 , where £ e is the number of accelerated electrons per 
proton. Naively, for the standard assumption that £ e ~ 1, one would expect \& ~ 10 -2 — 10 _1 
for typical values of e e ~ 0.1 and 1CT 3 < €b < 0.1. It is conceivable, however, that only a 
small fraction of the electrons participate in the acceleration process (i.e. £ e <C l). 3 

A comparison with observations shows, however, that \1/ 3> 1 (and £ e)max <S 1) for the 
vast majority of GRBs (Figs. 1 and 2). This implies that GRBs could arise from synchrotron 
emission in the external shock only if the efficiency of electron acceleration in relativistic 
collisionless shocks is very low (£ e <C 1). An external shock origin might still be possible if 
the radiation process responsible for the gamma rays is other than than synchrotron radiation 
(e.g., Wang et al. 2006). Alternatively, the prompt gamma-ray emission might arise from 
completely different mechanism, such as internal shocks (e.g. Ramirez- Ruiz & Lloyd-Ronning 
2002). 

It should be noted that afterglow observations already provide interesting constraints on 
the efficiency of electron acceleration (Eichler & Waxman 2005). Current observations imply 
that the characteristic energy of accelerated electrons is comparable to the proton post- 
shock temperature. They also imply that the efficiency £ e is similar for highly relativistic 
and sub-relativistic shocks and plausibly suggest that £ e ~ 1. However, even values of £ e as 
low as ~ m e /m p cannot be ruled out, since currently testable afterglow predictions remain 
unchanged for (E iso ,n cxt ) — >• {E iso ,n cxt )/£ e and (e e ,e B ) -> £e(e e ,es) for any £ e in the range 
m e / m p < £ e < 1 (Eichler & Waxman 2005). 

Estimates of the energy in the afterglow shock from late time radio observations when 
the flow is only mildly relativistic and starts to approach spherical symmetry (often called 
"radio calorimetry" ; Frail, Waxman & Kulkarni 2000; Berger, Kulkarni & Frail 2004; Oren 
et al. 2004; Frail et al. 2005; Granot et al. 2005) typically yield ~ 10 51 ' 5 erg assuming 
£ e = 1. However, as noted by Eichler & Waxman (2005), afterglow observations actually 
constrain ^ e Ek rather than E^. The true kinetic energies at late times are thus given by E\. ~ 
10 52 ' 5 (£ e /0.1) -1 erg. The initial energy content of the outflow could be even larger due to 
early radiative losses (i.e., during the prompt GRB and early afterglow stages). It is difficult 
to accurately account for the magnitude of such losses, as they depend on poorly known 
questions about postshock energy exchange between protons and electrons. Nevertheless, a 
lower limit on the radiated energy is given by E 1 = _E 7iiso (l — cos^o) ~ -E 7i i S o#o/2 (additional 
energy may be radiated outside the observed photon energy range, or during the early 



3 It should be pointed out that in principle even £ e > 1 is possible, especially near i?dec, due to pair 
enrichment of the ambient medium from pair production by gamma-ray photons that are scattered on the 
external medium (Ramirez-Ruiz et al. 2002; Meszaros et al. 2001; Kumar & Panaitescu 2004; Bcloborodov 
2002; Thompson & Madau 2000). 
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afterglow). Other possible channels of energy loss are the escape of accelerated non-thermal 
protons from the blast wave (high energy cosmic rays) or the production of high energy 
neutrinos via pion decay. Any added losses would inevitably lead to a further increase in 
the energy requirements. Therefore, very small values of £ e would imply very large energy 
contents. 

Another test of the simple external shock model is provided by a comparison of the 
correlation found by Firmani et al. (2006) between the isotropic equivalent luminosity, the 
burst duration, and the peak energy, with that predicted by Eq. 6. In the cosmo logical 
frame of the GRB this correlation reads L lso oc Ep 62±om TQ® 5 A9±om , where T0.45 is defined 
by Firmani et al. (2006) to be the time during which 45% of the counts above background 
are measured (which is expected to scale linearly with T GRB = T 90 ). This is in disagreement 
with Eq. 6, which, for a reasonably small scatter in e]/ 2 (e e /£ e ) 2 , gives L iso oc E^Tq RB . 

It is natural to hope that the values of e#, e e , p and £ e are universal, since they are 
determined by the microphysics of the collisionless shock. However, the wide distribution of ^ 
values seen in Fig. 2 suggests otherwise. That is, in the simplest version of the external shock 
model, a large scatter in ej/ 2 (e e /£ e ) 2 is required. The presence of a significant number of non 
shock-accelerated electrons in the external shock (£ e <C 1) appears to be more prominent for 
the sub-sample of bursts found to be spectroscopically associated with a supernova (Fig. 2), 
most of which have a smooth temporal profile (Kaneko et al. 2006). The low values of £ e do 
not, however, increase the total energy requirements to unreasonable values for these events 
as they have rather low values of E iso . Under this interpretation, a wide range of shock 
microphysical parameters may be the rule, rather than the exception. 

In conclusion, observations of the prompt emission in GRBs with known redshifts, which 
are becoming far more accessible in the Swift era, can provide an important diagnostic of the 
external shock model. Current observational constraints do not allow for efficient electron 
acceleration in the external shock, if its synchrotron emission produces the observed prompt 
gamma-ray emission. Although there is no a priori reason to suspect that £ e should be large, 
1 would dramatically increase the total kinetic energy budget. 

We are especially grateful to Stan Woosley and Eli Waxman for discussions. This work is 
supported by IAS and NASA through a Chandra Postdoctoral Fellowship award PF3-40028 
(ERR) and by the Department of Energy under contract DE-AC03-76SF00515 (JG). 
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Fig. 1. — * as a function of T 90 
Ghirlanda et al. (2004) and for all 4 GRBs with spectroscopically determined SNe (open symbols) from 
Kaneko et al. (2006). 
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Fig. 2. — Histogram of the inferred values of £ e , m ax = for the GRBs in Fig. 1. 



